Statefinder diagnostic for the modified polytropic Cardassian universe 
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We apply the statefinder diagnostic to the modified polytropic Cardassian universe in this work. 
We find that the statefinder diagnostic is quite effective to distinguish Cardassian models from a 
series of other cosmological models. The s — r plane is used to classify the modified polytropic 
Cardassian models into six cases. The evolutionary trajectories in the s — r plane for the cases with 
different n and j3 reveal different evolutionary properties of the universe. In addition, we combine 
the observational H(z) data, the cosmic microwave background data and the baryonic acoustic 
oscillation data to make a joint analysis. We find that Case 2 can be excluded at the 68.3% 
confidence level and any case is consistent with the observations at the 95.4% confidence level. 

PACS numbers: 98.80.Es,95.35.+d,98.80. Jk 



I. INTRODUCTION 



Recent type la Supernova observations, along with 
other observations such as cosmic microwave background 
(CMB) and galaxy power spectra, support the fact that 
thepresent expansion of our universe is accelerating 
P, 0. In order to explain the accelerated expansion of 
the universe, cosmologists have tried to explore many 
cosmological models. A negative pressure term called 
dark energy is always taken into account, such as the cos- 
mological constant model with equation of state lode = 
Pde/pde = — 1 where pde and pde arepressure and den- 
sity of the dark energy, respectively Q, the quiessence 
whose equation of state ujq is a constant between -1 and 
-1/3 [![, and the quintessence which is described in terms 
of a scalar field 4> 0, @ . Some other candidates are con- 
structed in different ways, such as the braneworld models 
which explain the acceleration through the fact that the 
general relativity is formulated in 5 dimensions instead 
of the usual 4 pj], and the Cardassian models which in- 
vestigate the acceleration of the universe by a modifica- 
tion to the Friedmann-Robertson- Walker (FRW) equa- 
tion 0]. In these cases, the dark energy component is 
not involved, but the accelerated expansion can still be 
obtained. The quantity w e ff is usually described as an 
effective equation of state for these models, and can be 
expressed by the Hubble parameter H and its derivatives 
with respect to redshift z jij. 

As so many cosmological models have been developed, 
a discrimination between these contenders becomes nec- 
essary. In order to achieve this aim, Sahni et al. proposed 
a new geometrical diagnostic named the statefinder pair 
{r, s}, where r is generated from the scalar factor a and 
its derivatives with respect to the cosmic time t, just as 
the Hubble parameter H and the deceleration parame- 
ter q, and s is a simple combination of r and q [9(. The 
statefinder pair has been used to discriminate a series of 
cosmological models, including the LCDM universe with 
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a cosmological constant A and a cold dark matter term 
(CDM), the Chaplygin gas, the holographic dark energy 
models, the quintessence, the braneworld models and so 
on. Clear differences for the evolutionary trajectories in 
the s — r plane have been found. 

In this paper, we apply the statefinder diagnostic to 
the modified polytropic Cardassian universe. The origi- 
nal Cardassian model based on two parameters O m o and 
n was first suggested by Freese & Lewis @. It was gen- 
erated from a modification to the Fricdmann equation. 
Such a universe is spatially flat and accelerating today. 
But it involves no dark energy term and is dominated 
by merely matter and radiation. These Cardassian mod- 
els predict the same distance-redshift relation as generic 
quintessence models, although they generate from com- 
pletely different physical principles. A generalized Car- 
dassian model-the modified polytropic Cardassian uni- 
verse can be obtained by introdu cing an additional pa- 
rameter f3 into this model [13, [HI, Q3 , which reduces to 
the original model if /3 = 1. The distance-redshift re- 
lation predictions of generalized Cardassian models can 
be very different from generic quintessence models, and 
can be differentiated with data from surveys of Type 
la Supernovae such as Super Nova/ Acceleration Probe 
(SNAP). In all, the modified polytropic Cardassian uni- 
verse can predict more fresh physical information than 
the original Cardassian. It is worthy of more detailed 
discussions. 

In this work, we successfully classify the modified poly- 
tropic Cardassian models into six cases by the statefinder 
diagnostic. The cases with different n and (3 correspond 
to different evolutionary trajectories in the s — r plane. 
The fact that LCDM corresponds to a fixed point (0, 
1) in the s — r plane plays a significant role for our 
classification. Also, it is very important to find where 
the evolutionary trajectories start and end. Another key 
standpoint is whether the evolutionary trajectory has a 
crossing with LCDM. We also study the relation between 
(n, (3) and the crossing redshift Zq, at which the modified 
polytropic Cardassian universe intersects with LCDM in 
the s — r plane. We find that the modified polytropic 
Cardassian models can be distinguished from other inde- 
pendent cosmological models by the statefinder diagnos- 
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tic. 

In addition, we use the observational H(z) data de- 
rived from ages of the passively evolving galaxies [HI , the 
newly measured value of the CMB shift parameter 1Z [l4[ 
and the „4-paramcter which describes the baryonic acous- 
tic oscillation (BAO) peak [l5[ to make a combinational 
constraint. We assume a prior of Hq = 72 ± 8 suggested 
by the Hubble Space Telescope (HST) Key Project [3. 
From the confidence regions, we find that Case 2 is not 
consistent with the observations at the 68.3% confidence 
level and all the six cases do not conflict with the obser- 
vations at the 95.4% confidence level. 

This paper is organized as follows: In Sec. 2, we briefly 
review the modified polytropic Cardassian universe. In 
Sec. 3, we introduce the statefinder pair {r, s}. In Sec. 4, 
we apply the statefinder diagnostic to various Cardassian 
models. In Sec. 5, we make a combinational constraint 
on the parameters of the modified polytropic Cardassian 
universe. In Sec. 6, the discussions and the conclusions 
are given. 



II. THE MODIFIED POLYTROPIC 
CARDASSIAN UNIVERSE 

Measurements of CMB suggest a flat geometry for our 
universe [l?], EH- If we consider a spatially flat FRW 
universe, the basic equation can be written as 
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where G is Newton's universal gravitation constant and p 
is the density of summation of both matter and vacuum 
energy. Freese & Lewis [|[ proposed a model called the 
Cardassian universe by adding a term on the right side 
of Eq.Q, 



H" = 



8ttG 



(2) 



where n is assumed to satisfy n < 2/3 and p m is always 
taken as a contribution of only the matter (in this paper 
we do not plan to consider radiation). If n = 0, it is 
identical to the cosmological constant universe. If B = 0, 
it reduces to the usual FRW equation, but with density 
of only matter. Thus, it is easy to get a new expression 
of H from Eq.©, 

H 2 = H*[n m0 (l + z) 3 + (1 - Q m0 )(l + z) 3n ], (3) 



by using 



Pm = Pmo(l + zf = ^mOPc(l + zf 



(4) 



where p m o is current value of p m and p c = 3Hq / 8t:G 
is the critical density of the universe. Clearly, this 
model predict the same distance-redshift relation as the 
quiessence with ujq = n— 1. But we notice that they have 
completely different essentials. The quiessence requires a 
dark energy component while the Cardassian does not. 



In the earlier time, the universe is dominated by the 
first term in Eq.© and the ordinary FRW behavior 
works. The additional term which consists of only matter 
gradually becomes a dominant driver afterwards. This 
transition was found to occur at z ~ 0(1) 0. Then the 
universe is caused to accelerate. The period of accelera- 
tion for this model is usually called the Cardassian era. 
The Cardassian model is attractive because the universe 
is flat and accelerating today but no vacuum energy is in- 
volved. And it has been demonstrated compatible with 
a series of observational tests, including the CMB data, 
the age of the universe, the structure formation and the 
cluster baryon fraction Q. 

The modified polytropic Cardassian universe can be 
obtained by introducing an additional parameter (3 into 
the above model [lO, EL LL2I , 
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H$[Sl m0 (l + zf + (1 - O m0 )/x(z)], (5) 
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(6) 

The two parameters (n,/3) are usually taken as n < 2/3 
and P > 0. If P = 1, the model reduces to the original 
one characterized by Eq. J3]) while if fx{z) = 1, the model 
just corresponds to LCDM. 

Similar to Wang et al., we also take fi m o = 0.3 as a 
prior in subsequent discussions. In the work of Wang 
et al., this model was compared with current supernova 
data and CMB data. It was proved that the existing 
data can be well fit for several chosen values of n and 
p. Also, the simulated data were constructed to make 
a discrimination between the modified polytropic Car- 
dassian universe and LCDM as well as the quintessence. 
Once fi m o is known with an accuracy of 10%, SNAP 
can determine the sign of the time dependence of dark 
energy density which provides a first discrimination be- 
tween various cosmological models fl2j . In this work, we 
use a geometrical tool-the statefinder diagnostic to make 
a classification and a discrimination about the modified 
polytropic Cardassian models. 



III. A BRIEF OVERVIEW OF THE 
STATEFINDER DIAGNOSTIC 

The Hubble parameter H = a/a and the deceleration 
parameter q = —a/aH 2 are two traditional geometrical 
diagnostics. They only depend on the scalar factor a and 
its derivatives with respect to t, i.e., a and a. Through 
a = 1/(1 + z), the deceleration parameter q can be ex- 
pressed as 



q(z) = ^(l + z)-l, 



(7) 



where H' is the derivative of H with respect to redshift 
z. The deceleration parameter q is a good choice to de- 
scribe the expansion state of our universe, but it is not 
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perfect enough to characterize the cosmological models 
uniquely. This shortage can be easily seen from the fact 
that many models may correspond to the same current 
value of q. And this difficulty can be overcome by an- 
other geometrical diagnostic-the statefmder pair {r, s}. 
This approach has been used to distinguish a series of 
cosmological models successfully. 

For a spatially flat universe, the statefmder r is defined 
as follows [i| 



the statefmder diagnostic in several literatures 0, 0| ■ In 
the s — r plane, a vertical line with s = n and r changing 
monotonically from 1 to 1 + 9n(n — l)/2 represents the 
evolutionary trajectory of the universe. 

For cases with ft ^ 1, we first pay our attention to 
the epoches in the far past (a — > 0) and the far future 
(a — > +oo). It is clear that z — > +oo represents the 
former case and z — > — 1 stands for the latter. From 
Eq. (ri~2"|) , it is easy to find the limit condition 



(8) 



where a is the third derivative of a with respect to t. s 
is just a combination of r and q, 
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3(9-1/2) 



(9) 



The statefmder pair was first introduced to analyze a flat 
universe with a cold matter and a dark energy term. For 
these contenders, r is given by 

9 3 Co 

r = 1 + -w D b(1 + w D e)^de - 2~H~^ DE ' ^ 

where J7rjE = Pde/Pc and wde is the derivative of wde 
with respect to t. And the other diagnostic, 



1 + uj - 1 ^ DE 

3 WDE-ff 



(11) 



From the two equations above it is easy to realize that 
LCDM corresponds to a fixed point (0, 1) in the s — r 
plane and the standard cold dark matter (SCDM) uni- 
verse with no dark energy term locates at (1, 1) forever. 
For this particularity, the current values of {r, s} pro- 
vide a considerable way to measure the distance from a 
specific model to LCDM. 

More generally, r and s can be given in terms of the 
Hubble parameter H and its first and second derivatives 
H 1 and H" with respect to redshift z, 

r(z) = l-2^(l + z) + ^ + (|^) 2 , (12) 



-2H'{l + z)/H + H"/H+{H'/Hf 
[) ~ 3(H'(l+z)/H-3/2) ■ 

Thus we can use the new tool to describe the evolutionary 
trajectories of the modified polytropic Cardassian uni- 
verse. 



IV. THE s-r PLANE FOR MODIFIED 
POLYTROPIC CARDASSIAN UNIVERSE 

For the original Cardassian model with ft — 1, it cor- 
responds to a compatible expression with the quicsscncc 
with luq = n — 1. The quicsscncc has been studied using 



lim r(z) = 1, 

z — >+oo 



(14) 



which means that the value of r in the far past is inde- 
pendent on n and f3. However, wc can not derive the 
similar properties for the value of s from Eq. (|13p . It is 
related to both of n and j3. For the limit of z —> — 1, from 
Eq.([12]), we get 



lim r(z) = 1 + -n(n — 1). 

2^ — 1 2 



And from Eq. lfTB")) . we have 



lim s(z) = n. 



(15) 



(16) 



The values of both s and r are independent on (3. 

Generally, some models have a crossing with LCDM in 
the s — r plane. By substituting the Hubble parameter, 
Eq.©, in the expressions of r and s, Eq. (fT2|) and Eq. (fl"3|) . 
we find that the crossing with LCDM happens at the 
redshift 



(3 + nP 



]3(l-»)/3 - 1. (17) 



And we notice that the crossing can happen only if the 
following inequality is satisfied, 



(fU - 1) 



1 - P + n/3 



> 0. 



(18) 



Due to the prior Sl m0 = 0.3 and /? > 0, — 1 > is 
naturally satisfied. Thus we may describe the condition 
in Eq. (fl"S|) cquivalently as 



1-/3 + 71/3 



< 0. 



(19) 



In order to understand the relation among n, (3 and 
zc clearly, we draw a f3 — zq plane for several fixed n in 
Figffl For n > 0, the crossing takes place at zc > for 
any ft, and zq changes little for larger values of ft. For 
n < 0, whether zq > or zq < depends on the values 
of ft. And 2c is nearly equal to -1 if ft is small enough. 
For the particular case which satisfies 1 — ft + nft = 0, the 
expected crossing happens at zc ^ oo. Clearly, n = 
and ft = 1/(1 — n) are the two critical conditions. 

Now we use Eq.(|19p to classify the modified polytropic 
Cardassian universe. We draw an n — ft plane in Fig [2] 
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FIG. 1: 0-Zc planes for n=-0.5, -1.0, -1.5 and -2.0 (the top 
panel) as well as n—0.2, 0.3, 0.4 and 0.5 (the bottom panel). 
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FIG. 2: n — (3 plane based on n and (3. The vertical solid 
line represents n = (Case 3), and the declined solid curve 
stands for (3 = 1/(1 - n) (Case 6). 



and divide it into different regions by the critical curves 
n = and (3 = 1/(1 - ri): 

■ Case 1 n< and < /3 < 1/(1 - ri); 

Case 2 n < and > 1/(1 - ri); 

Case 3 n = 0; ( , 

' Case 4 < n < 2/3 and /3 > 1/(1 — n); 1 J 

Case 5 < n < 2/3 and < /? < 1/(1 - ri); 
, Case 6 /3 = 1/(1 - n). 

Case 1 (n < and < (3 < 1/(1 - n)): n=-0.5, 
/3=0.1, 0.2, 0.3, 0.4 and 0.5 are taken for a qualitative 
analysis. We plot the s — r plane in FigJ3] Although 
the original Cardassian model with (3 = 1 is not involved 
in this case, we still draw this curve for a comparison. 
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FIG. 3: s — r plane for Case 1. The star (0, 1) corresponds 
to LCDM, the star (1, 1) SCDM and other dots the current 
values of the Statefinder pair, i.e., (so, ro) (also for subsequent 
cases). And the solid curves from bottom to top correspond 
to the evolutionary trajectories of n=-0.5, ,0=0.1, 0.2, 0.3, 
0.4 and 0.5 respectively. The dot-dash curves represent the 
critical case (3 = 2/3. For a comparison, the vertical dashed 
line stands for the case of (3 = 1. 



All the curves start at r = 1 and < s < 1 on the 
horizontal line, i.e., on the right of LCDM. After passing 
by LCDM, they arrive at their common end (—0.5,4.38) 
in far future. And the crossing with LCDM happens at 
z c =-0.999, -0.84, -0.50, -0.13 and 0.25 for 0=0.1, 0.2, 0.3, 

0. 4 and 0.5 respectively. Particularly, the fixed point (0, 
1) for LCDM just is the beginning point for the critical 
case of/3= 1/(1 -ri) = 2/3. 

Case 2 (n < and (3 > 1/(1 — n)): To be consistent 
with Case 1, we consider n = —0.5 again, but /3=0.8, 

1, 2, 3, 4 and 5 respectively The s — r plane is plotted 
in Fig[U The crossing with LCDM never occurs for this 
case. The vertical line s — n for the evolutionary trajec- 
tory of = 1 divides the whole plane into two parts. The 
left part with n < s < stands for 1/(1 — ri) < j3 < 1 and 
the right part with s < n corresponds to (3 > 1. And all 
the evolutionary trajectories commence at one point with 
r = 1. The same as Case 1, they arrive at (—0.5,4.38) 
in the end. 

Case 3 (n = 0): We plot the s — r plane for f3=0.3, 
0.7, 1, 2, 3 and 4 in Fig[5]. One point with r = 1 acts as 
the beginning point, and (0, 1) the common end. This 
diagram is divided into two segments by the LCDM fixed 
point, or equivalcntly the modified poly tropic Cardassian 
model with n = and [3 = 1. The segment with r > 1 
corresponds to (3 > 1 while the other segment with r < 1 
corresponds to (3 < 1. 

Case 4 (0 < n < 2/3 and f3 > 1/(1 - n)): The 
s — r plane for n = 0.2, (3=2, 3, 4 and 5 is plotted in 
FigJH] Clearly, all the curves start with r = 1 and s < 0, 
and then pass by LCDM after an arc route. All the lat- 
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FIG. 4: s — r plane for Case 2. The solid curves from right 
to left correspond to the evolutionary trajectories of n=-0.5, 
/3=0.8, 2, 3, 4 and 5 respectively. And the vertical dashed 
line stands for the case of j3 = 1. Same as Fig[3] the dot-dash 
curve represents the critical case /3 = 2/3. 
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FIG. 6: s — r plane for Case 4. The solid curves from right to 
left stand for the evolutionary trajectories of n — 0.2, /3=2, 
3, 4 and 5 respectively. The dot-dash curve represents the 
critical case j3 = 5/4. For a comparison, the vertical dashed 
line stands for the case of f3 — 1. 
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FIG. 5: s — r plane for Case 3. The solid curves from bottom 
to top stand for ,3=0.3, 0.7, 2, 3 and 4 respectively. The case 
of (3 = 1 just corresponds to the fixed point for LCDM, i.e., 
(0, 1). 



ter evolutionary parts of the trajectories nearly overlap 
with each other. They are insensitive to (3. The point 
(0.2,0.28) is the common end. The crossing redshifts 
are not too far from each other for different values of j3, 
i.e., zq ~ 0.3. It is interesting that the fixed point (0, 
1) for LCDM is just the beginning point for the critical 
case of /3 = 1/(1 — n) = 5/4. And the Cardassian mod- 
els for this case can satisfy the weak energy condition 
w = Px/{pm + Px) > ~ 1 although the effective equation 
of state satisfies ui c s = px/px < — 1- w c g < — 1 is con- 
sistent with many observations such as CMB and large 



FIG. 7: s — r plane for Case 5. The solid curves from right 
to left correspond to the evolutionary trajectories of n=0.2, 
/3=0.3, 0.6, 0.9 and 1.1 respectively. And the vertical dashed 
line stands for the case of j3 = 1. Same as Fig[6] the dot-dash 
curve represents the critical case /3 = 5/4. 

scale structure data (l9l I20I ] . 

Case 5 (0 < n < 2/3 and < /? < 1/(1 - n)): Same 
as Case 4, n = 0.2 is considered, but /3=0.3, 0.6, 0.9, 
1.0 and 1.1. We plot the s - r plane in PigEJ All the 
evolutionary trajectories start with r = 1 and s > 0. 
They arrive at (0.2, 0.28) in the end, and never pass by 
LCDM. The whole plane is divided into two parts by the 
vertical line s = n (corresponding to the case of (3 = 1). 
The left part satisfies 1 < (3 < 5/4 and (3 < 1 is satisfied 
for the right. 
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The Hubble parameter H{z) depends on the differen- 
tial age of the universe in this form 



FIG. 8: s — r plane for Case 6. The point (0, 1) stands for 
n = 0. Then n increases by a interval of 0.1 until 0.6 for 
cases on the right of (0, 1) and n decreases by a interval of 
0.1 until -0.5 for cases on the left of (0, 1). The other vertical 
dashed lines stand for the evolutionary trajectories of cases 
with (3 = 1 and different values of n. 



Case 6 (/3 = 1/(1 — n)): This is another critical case 
besides Case 3. n=-0.5, -0.4, -0.3, -0.2, -0.1, 0, 0.1, 
0.2, 0.3, 0.4, 0.5, and 0.6 are considered. For this case, 
crossing with LCDM happens at zq — > oo. We plot the 
s — r plane in FigO All the evolutionary trajectories 
start at LCDM in the past and end at (—1 + 3n/2, 1 + 
9n(n — l)/2) in the future. 

We have successfully given a qualitative analysis to 
the modified polytropic Cardassian universe from the 
statefinder viewpoint and it has been classified into six 
cases. A prior Q m o = 0.3 is adopted in this work. We also 
take some other values of f2 m o for comparison and notice 
that f2 m o is not a sensitive parameter for our analysis in 
this work. And it has been clear that the s — r plane is 
effective enough to discriminate the modified polytropic 
Cardassian models. The fact that LCDM corresponds 
to a fixed point (0, 1) in the s — r plane plays a signifi- 
cant role for our classification. Whether the evolutionary 
trajectory has a crossing with LCDM is our basic stand- 
point. Also, the beginning points and the ending points 
are both key factors for shapes of the evolutionary tra- 
jectories. 



V. DATA ANALYSIS FROM OBSERVATIONAL 
H{z) DATA, CMB AND BAO 

In order to understand the six cases above more clearly, 
we use the observational H(z) data derived from the 
passively evolving galaxies, the cosmic microwave back- 
ground (CMB) shift parameter 1Z and the ,4-parameter 
which describes the baryonic acoustic oscillation (BAO) 
peak to make a joint analysis. 



H(z) 



1 dz 



(21) 



1 + zdt 

which provides a direct measurement for H{z) through a 
determination of dz/dt. By using the differential ages of 
passively evolving gala xies determined from the Gemini 
Deep Deep Survey |2l| and archival data @, H, [H, [H[ , 
Simon et al. determined a set of observational H{z) data 
in the range < z < 1.8 and used them to constrain the 
dark energy potential and its redshift dependence [l3T |. 
Yi & Zhang first first used them to analyze the holo- 
graphic dark energy models and got a consistent result 
with others pt| . 

The model-independent shift parameter 1Z can be de- 
rived from CMB data. It is defined as fl4fl 



K = v 7 ^ 



ml) 



dz 



(22) 



E(zY 

where E(z) = H(z)/H and z r = 1089 is the redshift 
of recombination. From the three-year result of WMAP 
H3, Wang & Mukherjee estimated TZ = 1.70 ± 0.03 

Using a large spectroscopic sample of 46748 lumi- 
nous red galaxies covering 3816 square degrees out to 
z = 0.47 from the Sloan Digital Sky Survey, Eisen- 
stein et al. [l5[ successfully found the acoustic peaks in 
the CMB anisotropy power spectrum, described by the 
model- independent „4-paramctcr, 

'mol „L s f 1 ^] 2/3 , (23) 



where zi = 0.35 is the redshift at which the acoustic scale 
has been measured. Eisenstein et al. [l5| suggested the 
measured value of the ,4-parameter as A = 0.469 ±0.017. 

The best-fit parameters of the modified polytropic Car- 
dassian universe can be determined by minimizing 

v [H th ( Zi ) - H oh ( Zi )] 2 (rc- 1.70) 2 (.4- 0.469) 2 
^ a? + 0.03 2 



X 



0.017 2 

(24) 

where i?th( z i) is the theoretical Hubble parameter at z,-, 
H h(zi) is the observational Hubble parameter at Zi and 
Oi is the corresponding la error. We get the best-fit val- 
ues n = -1.85 and j3 = 0.23, with Xmin = !0-12. The 
best- fit results correspond to Case 1. And the cross- 
ing with LCDM in the s — r plane occurs at zc=0.32. 
The current values of the diagnostics arc so=-0.19 and 
ro=1.59. The confidence regions in the n — (3 plane are 
plotted in Fig[51 from which we find that Case 2 can 
be excluded at the 68.3% confidence level and all the six 
cases are consistent with the observational data at the 
95.4% confidence level. 



VI. DISCUSSIONS AND CONCLUSIONS 

The statefinder diagnostic is powerful to discriminate 
various cosmological models. Differences of the evolu- 
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FIG. 9: Confidence regions in the n — j3 plane for the joint 
analysis (the regions from inner to outer stand for confidence 
levels at 68.3%, 95.4% and 99.7% respectively). The vertical 
solid line represents n — (Case 3), and the declined solid 
curve stands for f) = 1/(1 — n) (Case 6). 



tionary trajectories in the s — r plane among a series 
of cosmological models have been found. For example, 
LCDM corresponds to a fixed point (0, 1) and the point 
(1, 1) represents SCDM. For the holographic dark energy 
model with c = 1 , the curve in the s — r plane com- 
mences at (2/3, 1) in the past and ends at LCDM in the 
future, with s monotonically decreasing from 2/3 to and 
r first decreasing from 1 to a minimum value and then 
rising to 1 [301 ] . Both the quintessence tracker models 
(with tracker potentials V = Vq /<f> a ) and the Chaplygin 
gas models have arc evolutionary trajectories, but in dif- 
ferent regions [J, [3l|, H3- The conditions — 1 < s < 
and r < 1 arc satisfied for the former cosmological model 
while < s < 1 and r > 1 for the latter. The evo- 
lutionary trajectories of the coupled quintessence models 
have more complicated evolutionary properties [33] . And 
quintessence models with other potentials were studied 
by Evans et al.[3l|, also along with a generalization of 
{r, s} to a non-flat universe. 

Also, the q — r plane has been widely used for discus- 
sion on the evolutionary property of the universe. For 
example, the point (0.5, 1) corresponds to SCDM in the 
q — r plane and the horizontal straight line from (0.5, 
1) to (-1, 1) stands for LCDM. The braneworld models 
have been studied too, includingDisappcaring Dark En- 
ergy (DDE) as the simplest case Q. LCDM separates the 
first braneworld model (named BRANE1 in [J] and the 
effective equation of state satisfies Ct> ff < —1) from the 
second branewor Id model (named BRANE2 in [| and the 
effective equation of state satisfies ui c g > — 1) and DDE 
models. DDE both begins and ends at SCDM, forming 
a loop. However, although both BRANE1 and BRANE2 
commence their evolutions at (0.5, 1) and end at (-1, 1), 
the diagnostic r satisfies r > 1 for BRANE1 while r < 1 



for BRANE2. 

As the s — r plane has been found robust for our clas- 
sification for the modified polytropic Cardassian models, 
wc do not intend to make use of the q — r plane. For 
the modified polytropic Cardassian universe, the evolu- 
tionary trajectories can be picked out easily with help 
of the s — r plane if f3 ^ 1. In fact, the original Car- 
dassian model with can not be discarded from the 
quiessence with ujq = n — 1 because they correspond to 
identical expressions for both r and s. Such a consistence 
has also been mentioned by Freese & Lewis, as well as 
how to distinguish the two models Q . 

Distinct differences in the s — r plane have been real- 
ized for the cases with different n and [3 for the modified 
polytropic Cardassian universe. For Case 1 and Case 4, 
the crossing with LCDM happens at some zc, while the 
same state never occurs for Case 2 and Case 5. And 
Case 3 and Case 6 just are two critical cases. Also, the 
beginning points and the ending points are tightly related 
to the shapes of the evolutionary trajectories. They se- 
riously depend on n and /3, especially n. As n and (3 
are found to be sensitive to the modified polytropic Car- 
dassian models, constraining the two parameters exactly 
becomes a valuable task. We use the observational H(z) 
data, the CMB data and the BAO data to make a com- 
binational constraint. We find that Case 2 can be ex- 
cluded at the 68.3% confidence level and all the six cases 
are consistent with the observational data at the 95.4% 
confidence level. Recent other constraints suggest some 
results far from consistent with each other. For example, 
in the work of Wang ct al., the choices of n = 0.2, j3 = 1; 
n = 0.2, (3 = 2 and n = 0.2, (3 = 3 are all consistent 
with SN la and CMB observations [13]. The first case 
corresponds to the quiessence with the equation of state 
cjq = —0.8, and this case is consistent with Case 5. The 
latter two cases with the crossing at zc=0.29 and 0.26 
respectively are in agreement with Case 4. And the cur- 
rent values of the diagnostics are so=0.13, ro=0.58 and 
s =0.16, r =0.45 for n = 0.2, /3=2 and n = 0.2, /3=3 
respectively. Nesseris & Perivolaropoulos [34[ suggested 
another result n = -23^, (3 = 0.025tj^? with a super- 
nova data set consisting of 194 SN la [35|, |3f|- Both Case 
1 and Case 2 may be included within the given error- 
bar. The best fitting values lie in Case 1 and the cross- 
ing with LCDM happens at zc=0.37. Meanwhile, the 
present quantities are so=-0.29 and ro=1.94. Evans et 
al. provided a best fitting result n = —0.94 and (3 = 0.06 
with a simulated data set [3l| . This result is in agreement 
with Case 1. The crossing happens in the far future, 
i.e., at 2c=-0-999. And the current values so=0.76 and 
ro=0.54 indicate a large distance from LCDM. As con- 
straints still remain weak, we expect for data with higher 
precision to provide more consistent fitting results in fu- 
ture. We also hope theses statefinder parameters can be 
determined more exactly and shed light on the nature of 
the cosmological models. 
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